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Ginzburg-Landau description of laminar-turbulent oblique band formation in 

transitional plane Couette flow 

Joran Rolland* Paul Manneville^ 
January 12, 2013 

abstract : Plane Couette flow, the flow between two parallel planes moving in opposite directions, is an example of 
wall-bounded flow experiencing a transition to turbulence with an ordered coexistence of turbulent and laminar domains 
in some range of Reynolds numbers [R s , Rt]. When the aspect-ratio is sufficiently large, this coexistence occurs in the form 
of alternately turbulent and laminar oblique bands. As R goes up trough the upper threshold Rt, the bands disappear 
progressively to leave room to a uniform regime of featureless turbulence. This continuous transition is studied here by 
means of under-resolved numerical simulations understood as a modelling approach adapted to the long time, large aspect- 
ratio limit. The state of the system is quantitatively characterised using standard observables (turbulent fraction and 
turbulence intensity inside the bands). A pair of complex order parameters is defined for the pattern which is further 
analysed within a standard Ginzburg-Landau formalism. Coefficients of the model turn out to be comparable to those 
experimentally determined for cylindrical Couette flow. 

1 Introduction 
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In their way to turbulence, wall-bounded shear flows display cohabiting turbulent and laminar regions. This striking 
phenomenon can even be statistically permanent and spatially organised, as for the flow between counter-rotating cylinders 
(cylindrical Couette flow, CCF) or counter-translating plates (plane Couette flow, PCF, Fig.[TJ top-left). Cohabitation then 
takes the form of alternately turbulent and laminar oblique bands. This peculiar pattern was first discovered by Coles and 
Van Atta in CCF (barber-pole or spiral turbulence) [I], the corresponding domain in the control parameter space being 
next charted by Andereck et al. [SJ. These experiments were restricted to the observation of a single spiral arm due to 
limited aspect-ratio (the ratio of the gap between the cylinders to the perimeter). 

Later on, Prigent et al. [3] performed studies at larger aspect-ratios, which allowed them to observe several intertwined 
spiral arms and to show that the oblique bands in plane Couette flow (Fig. [TJ bottom-left) were, qualitatively and quanti- 
tatively, the zero curvature limit of the spirals: Upon appropriate definition of a Reynolds number R based on the nominal 
shear rate, (i) these patterns bifurcate continuously at similar values of a well-defined upper threshold Rt above which 
turbulence is featureless, (ii) the spirals/bands are observed upon decreasing R down to comparable values of a lower 
stability threshold R g below which laminar flow eventually prevails, and (in) the streamwise and spanwise wavelengths are 
similar [4]. Figure [T] (top- right) recapitulates the experimental findings for PCF. 

Direct numerical simulations (DNS) of the Navier-Stokes equations for PCF were performed by Barkley & Tuckerman [5] 
- - -'who could obtain the band patterns in fully resolved, elongated but narrow, tilted domains. Their choice of boundary 
conditions however precluded the occurrence of patterns with defects or orientation changes inside the flow. This was not 
the case of the DNS by Duguet et al. [5] who recovered the experimental findings of Prigent et al. in fully resolved very 
large aspect ratio domains. Similarly, the spiral regime was numerically obtained by Meseguer et al. [7J and Dong [5] in 
CCF and the oblique band pattern in plane channel flow by Tsukahara et al. [9]. 

Up to now, there is no clear physical explanation for the formation of the spirals/bands from the featureless turbulent 
regime when R is decreased below Rt pUl b]. We however do have a consistent phenomenological description of the 
transition in CCF by Prigent et al. [3] in terms of two coupled Ginzburg-Landau equations with (strong) external noise 
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Figure 1: Top-left: Geometry of the plane Couette flow experiment, 2h is the gap, L x and L z ) the streamwise and spanwise 
dimensions, +U and — U the wall speeds. The Reynolds number is defined as R = Uh/v where v is the kinematic viscosity. 
Top-right: Experimental bifurcation diagram after Prigent [3, (a)]. Bottom: Picture of experimental banded turbulence in 
plane Couette flow apparatus with L x = 770h, L z = 340h at R = 358 (courtesy Prigent) and colour level representation 
of the perturbation energy averaged over y in our under-resolved ChannelFlow simulations with periodic boundary 
conditions at L x — 432, L z = 256, R = 290 and t = 18000. The two pictures are at roughly the same scale and similar 
Reynolds numbers after correction for the transitional range i?-shift due to numerical under- resolution [13 . 

added, introducing two complex amplitudes, one for each possible pattern orientation. Most of the coefficients introduced in 
these equations could be fitted against the experiments. In a similar vein, Barkley et al. [10] introduced the phase-averaged 
amplitude of the dominant Fourier mode of the turbulent mean flow modulation b] as an order parameter for the PCF 
transition. The emergence of the bands was then identified from the position of the peak in the probability distribution 
function (PDF) of this order parameter, shifting from zero in the featureless regime to a nonzero value in the banded 
regime. 

In the present article, we come back to the quantitative characterisation of the patterns in terms of order parameters. 
In contrast with [3[T0] we consider a configuration that does not freeze the orientation and allows for defective patterns. 
We keep the general noisy Ginzburg-Landau framework introduced in [3] for CCF and validate the approach in terms 
of amplitude equations at a quantitative level for PCF by means of numerical experiments. We take advantage of our 
previous work where the recourse to under-resolved DNS using Gibson's public domain code ChannelFlow [11] was 
introduced |12j . In |13) we brought evidence that this procedure could be viewed as a consistent systematic modelling 
strategy permitting simulations in wide domains during long time lapses at moderate numerical load. We indeed showed 
that all qualitative aspects of the transitional range are preserved at the recommended resolution (Fig. [T] bottom-right) 
and that, in the (slightly better) numerical conditions chosen here, the resolution lowering amounts to a 15-20% downward 
shift of [R g , Rt] from the experimental findings. This resolution reduction will allow us to accumulate statistics on moderate 
aspect ratio systems during very long times. We surmise that our results can be carried over to the realistic case of fully 
resolved simulations or experiments up to an appropriate adaptation of the Reynolds scale. We shall support this point of 
view briefly in £13.41 

We first recall the numerical procedure in §2.1| next we turn to the extraction of the turbulent fraction and the turbulence 
intensity ( §2.2p and to the definition of order parameters able to include information about the spatial organisation, 
£12.31 Results are then analysed in the successive subsections of £j3] devoted to the determination of the phenomenological 
parameters introduced by the Ginzburg-Landau formalism and accounting for the spanwise, streamwise and R dependence 
of the pattern. Section [4] summarises our findings. 
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Figure 2: Top: u x in an x, y plane, L x — 128 and in a z,y plane, L z = 180. Bottom: From left to right, u 2 in the 
V = —Vm plane, coarse-grained u 2 in the y < domain, resulting B/W discrimination, and W/G/B discrimination (see 
text). L x x L z = 128 x 180, R = 315. 

2 Simulations and data processing 
2.1 Numerical implementation 

The geometry of the experiment is described in Fig. [T] top- left. The Navier-Stokes equations are written in a reference frame 
where x, y, z are the streamwise, wall-normal, and spanwise directions respectively. Velocities are made dimensionless with 
U the absolute value of the speed at the boundaries y = ±h. Lengths are rescaled by h and time by h/U. The main control 
parameter is the Reynolds number R = Uh/v, where v is the kinematic viscosity but the flow regime also depends on the 
aspect ratios defined as T xz — L x , z /2h, where L x ^ z are the lateral streamwise and spanwise dimensions. In the numerics, 
h = 1 and the aspect ratios are T x _ z — L x ^ z /2. The base flow is independent of R : Vb = ye x . Written for the perturbation 
to the base flow u = v Vb , the Navier-Stokes equations read: 

dtUi + d 3 (uiUj) + yd x Ui + u y 8 i:X 

= - dip + R~ l d^Ui, 
d jUj = 0, 

with no slip boundary conditions at the plates, Ui{y — ±1) = 0, and periodic boundary conditions at distances L x and L z 
in the streamwise and spanwise directions, respectively. 

ChannelFlow [H] implements the Navier-Stokes equations using a standard pseudo-spectral scheme with Fourier 
transforms involving (N X ,N Z ) de-aliased modes in the streamwise and spanwise directions and N y Chebyshev polynomials 
in the wall- normal direction. As discussed in |13j . our numerical simulations are deliberately under-resolved: we use 
N y = 15, and N x<z /L XtZ — 8/3, which preserves all the qualitative features of the flow at a semi-quantitative level, just 
shifting the bifurcation thresholds down to R g = 275 ± 5 and i?t = 345 ± 5, to be compared with experimental or fully 
resolved numerical values, R g ~ 325 and Rt — 415 016]. 

In PCF, Prigent et al. experimentally found oblique turbulent bands with streamwise period X x ~ 110 and variable 
spanwise period A 2 from 85 around i? g to 45 close to Rt- The sizes of our numerical domains range from L z = 24 to 
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192 and from L x = 80 to 170. Our domains hence remain rather small since they can contain one to three such spanwisc 
wavelengths but they are much larger than the minimal flow unit |14) of size £ x ss 6 and £ z sa 4, below which turbulence 
cannot self-sustain. They are also much longer in the streamwise direction than the tilted domains considered by Barkley 
et al. 5, 10 but remain smaller than the largest domains considered by Duguet et al. [6] or in our preliminary studies |13] 
which went up to L x = 800 and L z — 356 but at a much lower resolution, or the latest experiments by Prigent et al. with 
L x = 770 and L z = 340 0. 

2.2 Local averaging and related quantities 

The square of the perturbation velocity u 2 is a good indicator of the local state of the flow. Figure [2] (top) displays colour 
level representations of that quantity in typical wall-normal planes, streamwise (x, y) with height 2 and length L x = 128, 
and spanwise (z,y) with height 2 and width L z = 180, for R = 315. The pattern seen from above in the (x,z) plane at a 
given wall-normal coordinate y = —0.57 = — y m is displayed in Fig. [2] (bottom, left), the other panels represent the same 
image after additional post-treatment to be discussed below. The value y m = 0.57 roughly corresponds to the place where 
u 2 is statistically the largest in the range of Reynolds numbers of interest (see Fig. 5 in |13)V 

The simplified representations shown in the centre and right panels of Fig. [2] rest on the coarse-graining of the u 2 field 
introduced in [12]. This procedure directly stems from the general organisation of the flow in the band regime already 
identified in previous studies p][5] and clearly visible in the side and front views of the flow in Fig. [2] (top). These 
pictures suggest to average over the upper layer of the flow (y > 0) and its lower layer (y < 0) separately. Typical 
experimental observations [2j|3][l5], film and pictures, yield an information integrated over the whole gap, which motivates 
us to compute comparable quantities. As shown in Fig. [3j the computational domain is divided in small stacked boxes of 
size l x xl y xl z = 2x1x2. This size is slightly smaller than (but related to) that of the minimal flow unit. The width l z = 2 
approximately corresponds to the spanwise size of a turbulent streak. By contrast, l x = 2 is much smaller than the typical 
length of a turbulent streak, I ~ 40, so that the turbulent intensity variations along a streak can be captured. Quantity u 2 , 
henceforth called 'energy' by a small abuse of language, is then averaged in each of these cells and a threshold c is chosen 
according to which it is laminar or turbulent. The turbulent fraction / is then the proportion of turbulent cells, and the 
turbulent energy et is the energy conditionally averaged in space over the turbulent zone. Conditional averaging of any 
field can easily be performed in the same way. The reduction procedure is expected to depend on the value of c. As seen 
in Fig. |4] which displays the profile of the coarse-grained energy through the band pattern, the locally turbulent flow has 
typical energy higher than 0.1 and locally laminar flow less than 0.05. The computation of the time-averagecfj] turbulent 
fraction F and the time-averaged turbulent energy E t for values of c ranging from 0.005 to 0.13 did not pointed to an 
optimal value for c, as expected from a flow displaying a smooth modulation of turbulence, and c = 0.025 was eventually 
chosen with little consequence on the quantitative information drawn from the procedure. 

A typical example of this thresholding is given in Fig. [2j bottom line: from a realisation of the flow at y = —y nl (left) we 
compute the coarse-grained energy for y < (centre-left) and apply the criterion to obtain a black-and-white (B=laminar, 
W=turbulent) representation of the flow, still for y < (centre-right). Distinguishing the y > layer from the y < layers 
allows a refined representation of the flow as shown in the bottom-right panel of Fig. [2] which displays the turbulent and 
laminar areas using a black/gray/white code: 'black' represents laminar cells of top of each other, 'white' turbulent cells on 
top of each other, 'light grey' y > turbulent cells on top of y < laminar cells, and 'dark grey' y < turbulent cells of top 
of y > laminar cells [T2"] . As already seen in the top panels, the streamwise direction going from left to right, turbulence 
is to the right of the band for y > and to its left for y < 0, in agreement with previous findings Q][S]. This fact could be 
used to compute properties at the edge of the bands, for instance velocity or energy profiles. A quantitative comparison to 
results of Barkley and Tuckerman [5j b] has not been attempted since the differences in geometry and resolution shift the 
Reynolds number correspondence. 

The procedure has been implemented on-line to allow the computation of time series of the turbulent quantities. Since 
these quantities fluctuate, we compute their time-averages E, E t , and F as 



where To is introduced to take into account the transient necessary for the flow to reach its permanent regime, and T is 
taken sufficiently large (typically, over 5000) to keep the relative fluctuations of E(T) within 0.5 %. The cut-off c being 

1 On general grounds, lower case letters will denote instantaneous values and upper case letters the corresponding time averages. 
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Figure 3: Sketch of the averaging boxes from the side and from above. 




Figure 4: Coarse-grained u 2 -profile projected along the direction of the turbulent band (L x x L z 
y<0). 



128 x 64, R = 315, 



appropriately chosen, the energy content of the laminar part is negligible so that we have e ~ / x et , which means that 
the average energy of the flow is positively correlated to the changes of turbulence intensity in the bands measured by 
et, as well as to the fractional area / occupied by the bands. On the other hand, quantities et and / do not show much 
correlation. This can be seen in figure [5] in which normalised quantities, e = (e — E)/((e 2 ) — E 2 ) 1 / 2 , etc., are displayed. 
Computation of the correlation of e and /, as well as e and e t yields 0.5 ± 0.1, on average over all experiments, whereas e t 
and / are not correlated, yielding ± 0.1. Owing to small relative fluctuations, the relation e ~ / x et implies a similar 
relation, E ~ F x E t , for the averaged quantities. 

It turns out that e(i), e t (t), and f(t) are little affected by the orientation fluctuations: even when the pattern presents 
defects, the surface occupied by turbulence and the turbulence intensity in the bands remains essentially unchanged. This 
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Figure 5: Time series of the normalised average energy e, turbulent energy et and the turbulent fraction / in a typical 
numerical experiment for L x x L z = 110 x 48 and R = 315. 



allows us to perform averaging regardless of the orientation, but E, Et, and F remain sensitive to the value of the bands' 
wavelength imposed by the periodic boundary conditions fixing the in-plane dimensions L x>x , as discussed below. 

2.3 Order parameter 

2.3.1 Conceptual framework and operational definitions 

In the theory of phase transitions, an order parameter is an observable which, at the thermodynamic limit (permanent state 
at infinite size), is zero in the non-bifurcated state, here the featureless turbulent regime, and non-zero in the bifurcated 
state, here measuring the amount of coexisting laminar and turbulent domains. The turbulent fraction F (introduced 
in |15] at a time when the spatially organised character of the banded regime was not yet recognised) or rather the laminar 
fraction 1 — F, partially fulfils this condition but remains of limited value since it does not account for the space periodicity 
of the pattern explicitly, which is what we want to overcome, inspired by previous work [3lll0j. In pattern- forming systems, 
the bifurcation is generally characterised by the amplitude of the relevant bifurcating mode and, especially in extended 
systems, by the amplitudes of the modes entering the Fourier decomposition of the structure that develops from the 
instability mechanism. When fitting the pattern-forming problem into the phase transition formalism, these amplitudes 
are the natural order parameters. 

Figure [6] illustrates the result of a Fourier analysis of patterns with three bands fitting a domain of size L x x L z — 128 x 
180. Symmetries in the spectrum allow us to consider wave numbers such that < n x < N x /2, —N z /2 + 1 < n z < N z /2^ 

2 Strictly speaking only < n x < N x /3 and —N z /3 + 1 < n z < N z /3 since N x ,z are the numbers of de-aliased modes so that, the 3/2-mlc 
being used, the number of modes truly involved in the dynamics is 2N XyZ /3 and the corresponding bounds (2N X:Z /3) /2 = N X}Z /3. This proviso 
is however not essential since we are only interested in centre of the spectrum with n x<z small. 
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Figure 6: Color plots of patterns and natural logarithm of the corresponding spectra averaged over y as explained in the 
text for a well-formed pattern (left) and for a pattern with defects (right). In the spectra, u x (0, y, 0) is set to zero, which 
yields a black dot in the spectra; normalisation by the maximum value makes its position appear the peak as a white dot. 
L x x L z = 128 x 180, R = 315. 

The figure displays (x, z)-plots of u 2 at y — ~y m (left) and corresponding spectra averaged over the wall-normal direction 
(right). The top panels correspond to an ideally formed pattern and the bottom panels to a defective one. For both 
flows, the wave numbers corresponding to the peak are n x — I and \n z \ — 3. The spectra are zoomed on the smallest 
wave numbers so that modulations at the scale of the streaks are outside the reframed graphs. When the pattern is well 
formed, a single mode corresponding to the fundamental of the modulation clearly emerges, about two orders of magnitude 
larger than the other modes. These background modes account for small irregularities at a given time and not to steady 
anharmonic corrections to a basically sinusoidal profile: the average ratios m„^3/m„ i= 3 are at most 0.1 and the harmonics 
have no definite phase relation with the fundamental, corroborating the observation by Barkley and Tuckerman that the 
modulation is quasi-sinusoidal b]. In the defective case (figure El bottom), two peaks emerge, corresponding to the two 
orientations. Their amplitude is smaller, and other harmonics have non negligible amplitudes, accounting for the spatial 
modulations of the pattern. Envelopes can be defined, one for each orientation, obtained by standard demodulation. 

The picture shown corresponds to a case with three bands, showing that there is enough room for a grain boundary. 
For smaller systems with one or two bands, defects correspond to the coexistence of laminar and turbulent regions without 
conspicuous organisation. In fact, the pattern can be observed only when the domain is above some minimal size L™™. 
Our simulations suggest £™ m ~ 24 and L™ ln ~ 70 (a precise determination of the minimal size is still under study). This 
is much smaller than in experiments because periodic boundary conditions tend to stabilise the pattern: only a tendency 
to form oblique turbulent patches was observed in laboratory experiments with L x x L z = 280 x 72 [15j . where the ideal 
simple shear flow was achieved with sufficient accuracy only in the centre of the set-up due to lateral boundary effects. 

Following Prigent et al [21 c], both orientations being equivalent, we expect that the pattern can be characterised by 
two complex quantities A± : 



where A± £ C describe slow modulations at scales much larger than X xz = 2ir/k x z , the 'optimal' streamwise and spanwise 
wavelengths. Variables x and z denote the corresponding space coordinates. Despite the highly fluctuating nature of the 
turbulent flow, the pattern being time- independent, there is just a possible slow evolution at an effective time t linked to 
wavelength selection and defect dynamics. The modulus of A± gives the amplitude of the turbulent intensity modulation, 
and the phase fixes the absolute position of the pattern in the domain. 

Near the threshold Rt, introducing e = (i?t — R)/Rt, A± are guessed to fulfil Ginzburg-Landau-Langevin equations in 
the form [3J: 




(1) 



± 



T dtA± 



-9i\A±\ 2 A± -9i\A T \ 



2 A± + a(± 



(2) 
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Figure 7: Top: Variation with time of the phase of the main Fourier component of u x for L x = 110, L z = 48, R = 290. 
Bottom: growth of the variance of the ensemble-averaged fluctuations as a function of time (T e = 130, N e = 1000), with 
linear fit. 

where the a£± are additive noise terms expressing the local fluctuations caused by intense small scale turbulence, a being 
the strength of the physical noise. Though this noise is both more intense than thermal fluctuations (see [TB] and references 
therein) and much more correlated since the featureless turbulent state is not without structure |18j . terms £± are tacitly 
taken as independent normalised delta-correlated Gaussian white noise processes ((C±(*)C±(*')) — ${t — *'))• 

Periodic boundary conditions determine accessible wavelengths in a given domain: X x _ z = L x _ z /n x ^ z , where the integers 
n IjZ are the wave numbers. In the computational domains considered here, with L X}Z not so large, it turns out that states 
with n x — 1 or 2 and n z = ±1 up to n 2 = ±3 can be observed, depending on the precise value of L x and L z . When the 
wavenumbers are small enough, the partial differential model ([2]) can be reduced to a set of ordinary differential equations 
for scalar complex amplitudes, and when there is no wavelength competition but only an orientation competition playing 
with the ±, just by two amplitudes ^4± jT i x ,n z corresponding to a specific pair of wavenumbers (n x ,±n z ). These amplitudes 
are then governed by: 

- (ffl|^±,n.,n,| 2 + S2|^,n.,n.| 2 )^±,n. > n.+aC± (3) 

with e nx ,n z = e — ^ x ^k^. — £ z 5k z , 5k x . z = k XiZ — k x z , and k XiZ — 2im X:Z / L X:Z , so that the dependence of the pattern on the 
value of the wavevectors can be studied by changing the size of the domain. 

When a single wavelength and a single orientation are selected, a single complex amplitude can serve to characterise the 
corresponding pattern. This was precisely the case considered by Barkley & Tuckerman [5] who defined the order parameter 
from a single Fourier amplitude by sampling its probability distribution function (PDF) and averaging over its phase [10) . 
So doing, they were able to detect the bifurcation to the band regime from the change in the PDF as i?t was crossed. In 
our simulations a single orientation is selected only deep enough in the band regime, i.e. sufficiently below R t but above 
R g . The pattern is then well installed and its orientation remains fixed but its lateral position in the domain can fluctuate, 
which strictly corresponds to the phase fluctuations alluded to above. When this is the case, symmetry considerations 
underlying ([3]) imply that the phase is dynamically neutral, hence constant in a deterministic context, while it is expected 
to evolve as a random walk in a noisy context [10]. Figure [7] shows that this is indeed the case. The top panel illustrates 
the variations of the phase of the Fourier amplitude u x of the streamwise velocity component at y — —y m for R = 290 in 
a domain of size L x x L z — 110 x 48. The expected property is illustrated in the bottom panel which displays the linear 
growth of the variance of the phase fluctuations as a function of time after appropriate ensemble averaging: From the initial 
time series we define an ensemble of N c successive sub-series of duration T c as: 

0i(t) = 0(t + (i-l)T e )-0((i-l)Te) , te [0,T e ] 
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for i = 1 . . . iV c . We next define the ensemble average: 

(#) = ^-&W- te[o,r e ]. 

z— 1 

which always remains of order 10~ 3 , while the variance: 

1 We 

4(*) = ]^E^(*)-W(*)) 2 ' *6[0,T e ] 

i= 1 

is indeed seen to grow linearly with time (Fig. [7]). 

When the wavelength and/or the orientation can fluctuate, as is now the case of interest, the practical definition of 
an order parameter is less straightforward since a single complex amplitude is not enough. Here, we forget about the 
information contained in the phase of the relevant complex amplitudes and focus on their modulus. We then define the 
instantaneous order parameter m n;Ej „ T (t) as the modulus at time t of the fundamental Fourier mode (n x ,n z ) accounting 
for the pattern as featured by the streamwise velocity field u x averaged along the wall-normal direction: 

(\ f 1 \ 1/2 
m n x ,nAt) = ( 2 J \u>x(n x ,y,n z ,t)\ dy) (4) 

but equivalent results are obtained from the other velocity components, with or without wall-normal averaging. 



2.3.2 Typical experiments and the order-parameter time-averaging issue 

Owing to the linear stability of the laminar flow, turbulence has to be triggered by finite amplitude disturbances. A 
typical experiment consists of creating a random initial condition and evolving it at a Reynolds number for which uniform 
turbulence is expected, here R = 450 (3> i?t — 345 at the resolution chosen in the present work). That state is next 
used as an initial condition for a simulation at the targeted value of R for which the pattern of interest is expected, hence 
R g < R < Rt- Such experiments were named quench in [15,17 . Variations of turbulent quantities e(t), et(t), f(t) and of 
the order parameters m(t) at the beginning of a typical experiment are shown in Fig. |8l The stabilisation of the featureless 
regime at R — 450 is clearly visible with e ~ et, //10 ~ 0.1 (left), and m ~ (right). The subsequent quench at t = 500, 
R = 315 is seen to produce some undershoot of e, et and /, while m grows slowly, which corresponds to the formation 
of bands. After a short period of exponential growth, the order parameters saturate as shown in Fig. [9] for a series of 6 
independent runs in the same conditions where a band is expected, pointing out the selection of the orientation, with one of 
the order parameters larger than the other by typically one to two orders of magnitude. The simulation is continued during 
at least 5000 time units in order to ensure good convergence of the time averages E, Et, F and M. The same procedure is 
repeated for all the values of L x , L z , and R considered, except in £13.41 where an adiabatic procedure is adopted to vary R. 

Like the turbulent quantities e, et and /, order parameters m„ xi „ z fluctuate in time but, since orientation changes are 
now of interest, care is required when computing their averages. Figure [TUl displays a typical example of long-lasting time 
series of m 2 for L x = 128, L z = 84, and R = 315, which produces patterns with n x — 1 and \n z \ — 1 and 2, so that 
modes n z — ±1 and ±2 dominate in turn. As long as the instantaneous state of the system is close to ideal, m fluctuates 
around a specific mean value which depend only on \n z \ as expected from symmetry considerations. Defects may appear 
and disappear, involving several modes with similar amplitudes. For the data in Fig. 1101 L z — 84 lies in a range L z G]80, 96 [ 
where the competition between different values of \n z \ is particularly intense (see below). When it is the case, a proper 
definition of order parameters implies conditional averaging over periods during which the pattern is well formed with the 
chosen value of \n z \. For example, in Fig. 1 1 01 \n z \ — 1 is present during about 3/4 of the time window and \n z \ = 2 less 
than 1/4 of it. Since m±i > m±i when the corresponding modes dominate the pattern, one gets Mi > Mi, but it would 
be meaningless to make a blend of the two and define a single order parameter for the system. A detailed study of this 
special case is deferred to [19] . 

However, outside cases of strong wavelength competition, a single value of \n z \ is selected, which makes things somewhat 
easier and allows us to simplify the notation: m nx} ± n ^ i— > m±. An example is displayed in Fig. [Tl] for L z = 32 where only 
\n z \ — 1 shows up. Averaging can then be performed from two-dimensional probability distribution functions (PDF) 
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Figure 8: Time series of turbulent quantities (left) and m (right) during the initial stage of a typical experiment. Here, for 
L x x L z = 110 x 48, R is initially set at 450 and switched to 315 at t = 500. 




Figure 9: Time series of m 2 for six different runs at R — 315, with L x x L z = 128 x 48, starting from uniform turbulence 
(time is reset upon quenching). + and — refer to n z = +1 and n z = —1. Mode +1 is selected in runs 2, 3, 4, and mode — 1 
in runs 1, 5, 6. (colour online) 
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Figure 10: Time series of m 2 for L x x L z = 128 x 84 at R = 315. Well-formed patterns with n z = +1 appear for t € [43, 51], 
n z = —1 for t € [22,37] and t > 52, n z = +2 for t G [3, 11], and n z = —2 for t G [14, 16]. Defective patterns are obtained 
for t G [18,22] or [38,43]. (All times to be multiplied by 10 3 .) 




Figure 11: Time serie of m 2 (t) for L x x L z = 110 x 32, R = 330 



x 1 0" 3 x 1 0" 3 




m m m m 



Figure 12: PDF of m for L x x L z = 110 x 32, R = 290 (left), R = 330 (centre-left), i? = 330, symmetrised (centre-right), 
R = 337, symmetrised (right). 
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II c (m + , m_), where subscript 'e' means 'empirical'!! Far away from R t the orientation does not fluctuates and the pattern 
remains without defects, which yields a one-hump PDF such as the one in Fig. [12] (left) for R — 290, but closer to Rt the 
orientation fluctuates and defects are present. Two humps are then obtained as in Fig. 1121 (centre-left) which derives from 
the time series for R = 330 shown in Fig. [TT] Due to the finite length of the time series, the PDF is not symmetrical with 
respect to the diagonal but since, for symmetry reasons, the two orientations should be present with the same weight, one 
may improve the statistics by constructing Yl CtS (m+, m_) = | (I1 c (to + , m_) + 1I g (to_, rn+)), where the additional subscript 
's' means 'symmetrised', which is done in Fig. [12] (centre-right). Averages can then be extracted from the 'symmetrised' 
PDF, which works fine as long as the orientation fluctuates but neither the wave numbers n x nor \n z \. We thus define: 

M = 2 m'n e>s (m', m") dm'dm" . 

J m" <m' 

The right panel in Fig. 12 displays the (symmetrised) PDF for R = 337, when re-entrant featureless turbulence intermittently 
bursts in, which manifests itself as a secondary hump close to the origin, see below £|3.4i and especially the discussion related 

to Fig. [m 



3 Results 

3.1 Theoretical expectations 

The statistically steady states (permanent regimes) obtained in the DNS and characterised by the time-averaged empirical 
order parameters to„ x „, defined through Q can then be compared to the equilibrium states predicted by model ([3]), the 
deterministic part of which can be written as deriving from the potential: 

V = -i?(|A + | 2 + |A_| 2 ) + i 9l (|A + | 4 + |A_| 4 ) 

+ \g 2 \A + \ 2 \A-\\ (5) 

where e is a short hand notation for e nxi „ 2 , computed from the values of n x and ±n z relevant to the pattern of interest, 
again with a single pair of modes present in the system. 

Assuming Gaussian noises of strength a, the theoretical expression of the PDF reads [20 : 

n t (m + , m_) = Z^ 1 m + m- cxp(— 2V/a 2 ) , 

where subscript 't' means 'theoretic' and 



OO p oo 



Z= / / 77i + m_ exp(— 2V/a )dm_)_dm_ 
Jo Jo 

is a normalisation factor called the partition function in statistical physics. For values of e that are not too small, the most 
probable values m^. corresponding to the maxima of n t give a good estimate of expected mean values (m±) (mean-field 
approximation). They are given by the solutions to: 

= — em ± + (?iTO ± + gim±m^ — a /2 . 

At lowest non-trivial order in a 2 , we have: 

m± oc — , (6) 

which corresponds to the trivial solution of the deterministic problem, just shifted by the effects of noise. The non-trivial 
solutions read: 

< = m- = \M<h+92) (7) 



3 In contrast with what was defined by Barkley et al. |10| who chose to scale out the pre-exponential factor, having dP(a) = adap(a), where 
a is the modulus of the dominant Fourier mode, corresponding to one of our m±, we have here dil(m+,m_) = dm+dm_ n(m+,m_), as a 
consequence: 11(0, m_) = il(m-|_,0) = 0. 
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Figure 13: Perturbation energy E (left), turbulent fraction F (center) and turbulent energy E t (right) as functions of L z 
for L x = 128 and R = 315. 



and 

m° ± = y/i/^, m° = = (8) 

V 2e (.92 - 9i)/9i 

For e < 0, solution ([6]) is stable and the other solutions do not exist. For e > 0, solution ([6]) is unstable and the symmetric 
solution ([7]) is a saddle point since, in order to get a stripe pattern we assume <?2 > 9i] otherwise a stable rhombic pattern 
would be obtained but is observed neither in the experiments nor the numerical simulations. This solution lies on the 
boundary of the attraction basins of solutions (|SJ), which exist for e > and are stable. They represent the amplitude of the 
turbulence modulation for R < R t . The mean amplitude of the installed mode varies as Ve as is typical of a supercritical 
bifurcation. The other mode, expected to be zero in the deterministic case, is present with small amplitude due to noise. 
Noise is also responsible for a switch from the '±' situation to the '=f' one when fluctuations make the system leave the 
well corresponding to an installed '+' mode to reach the other one where the '— ' mode is installed and vice versa, going 
through the potential barrier at the saddle solution (J7J . The asymptotic expressions above agree with the values computed 
from the PDFs obtained by direct simulations of model (|3]) and will be plotted together with our results in Figures QjJ] 
and 1211 When e is very small, fluctuations around the most probable values have to be taken into account. The mean field 
approximation is no longer valid and a behaviour in the form (m±) oc \e\P is expected, where /3 is the critical exponent 
describing the variations of the order parameter with the control parameter in the theory of phase transitions. We shall 
restrict to the mean-field approximation as a first guess since the nature and extent of this specific regime, called critical in 
statistical physics, are not yet clearly characterised in the present case (see [16] and references therein for examples where 
fluctuations have thermal origin). 

The deterministic part of model ([3} is invariant against phase changes of the complex amplitudes A± = m± exp(i(f>±), 
implying that the <f>± are dynamically neutral. They are indeed governed by: 

T di<j)± = 1m [exp(-i<p±)a(±] /m± , (9) 

i.e. a stochastic process, the strength of which depends on the instantaneous value of m±. In fact, the right hand side of 
(|9|) is another random Gaussian process Q(t)a/m{t) with zero mean and variance a 2 S(t — t')/((m} 2 + a 2 ^) where <r^ is the 
variance of m(t), which can be checked numerically using model ([3]). Results in Fig. [7] above can be quantitatively rendered 
by taking a/r = 4 x 10~ 4 . 

Coherence lengths £ a and £ z in control how strictly the wavevectors k x and k z are bound to their optimal values 
k% and k c z . The anisotropy of the base flow leads to expect different values for £ x and £ z . For PCF, experimental data [3] 
suggests that \ x and therefore k x do not depend on the Reynolds number, whereas X z decreases with R. Prigent et al. also 
report a decrease of the effective value of £ 2 as R is increased but the experiment did not give access to £ x . In the following, 
we determine most of coefficients in model (j3]) from the dependence of E, E t , F and M on k z , k x , by varying L x , L z and 
R using the quench protocol explained above. The dependence upon the Reynolds number analysed next is obtained from 
simulations in which R is varied adiabatically. 
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Figure 14: Order parameter squared M as a function of k z (left) and of (k z — k c z ) 2 (right) for L x = 128 and R = 315. 



3.2 Dependence on k z 

We fix R = 315, in the middle of the range where bands are expected at the resolution that we consider [13], and L x = 128 
so that a single streamwise period is obtained (n x = 1, k x = 2ir/L x ). We take values of L z ranging from 24 to 192. Taking 
the number \n z \ of spanwise periods into account, we have k z — 2-n\n z \/ L z . Figure [13j displays E, E t , and F as functions 
of L z , showing that \n z \ increases with L z : one band for 24 < L z < 80, two bands for 96 < L z < 144 and three bands for 
156 < L z < 192. In these ranges, \n z \ stays fixed during the simulation. In contrast, patterns with \n z \ = 1 and \n z \ = 2 
alternate in time for 80 < L z < 96, here for L z = 84 (see figure [TU]) and L z = 90. This special case is studied more 
thoroughly in [15]. A similar competition between \n z \ — 2 and \n z \ = 3 is expected to occur for 144 < L z < 150. Taken 
together, the results in Fig. [13] illustrate confinement effects when L z is small. Turbulence is featureless for L z < 24 and 
the turbulent fraction F (central panel) rapidly decreases from 1 down to ~ 0.63 which therefore represents some kind of 
optimum at R — 315. 

The results also suggest to check cases with n z > 1 against case n z = 1. Figure [T4l (top) displays M 2 as a function of 
k z = 2ttti z / L z and n z — 1 as a full line. Data obtained with two and three bands are also shown as dashed and dash-dotted 
lines, respectively. For them no points at large wavevectors are obtained because the corresponding patterns are not stable 
enough to be observed. The parabolic shape expected from the theory ( i)3.ip is reasonably well reproduced by the data. The 
maximum is reached at k\ ~ 0.16, that is A^ = 2ir/k c z ~ 39, as determined from a fit against a parabola. The so-obtained 
value of k z can next be used to determine £, 2 from the slope of a linear fit of M 2 against (k z — kf) 2 . The result is displayed 
in Figure Q3] (bottom) where data corresponding to one band are shown with '+' signs. From it one derives £, 2 /gi = 0.1. In 
turn, the constant term in the fit is a compound accounting for the dependence on R and k x , namely (e — £, 2 6k 2 )/gi = 0.002. 
Data for two and three bands, respectively shown with 'o' and 'x' symbols, are seen to be consistent with these estimates. 
Here a single value of R has been considered. In the CCF case, Prigent et at found for £, 2 /gi values of the same order of 
magnitude, decreasing with R from 0.5 to 0.1 [3J c]. 



3.3 Dependence on k x 

The dependence of the pattern's characteristics on k x is studied for R = 315, L z — 48, and L x g [80, 170]. In this range, 
only n x = 1 is obtained, except for L x = 170 where n x = 2 can also be observed. Figure [15] shows that, as a function of 
k x (top), M 2 displays a maximum at k x = 0.058, hence = 110, whereas fitting M 2 against (k x — k%) 2 (bottom) yields 
H/gi = 2.7. The same study at R = 330 (closer to R t = 345) gives A^ = 110 and = 3.9, while at R = 290 (closer 

to R g = 275) we get A£ = 125 and £, x /gi — 2.2, which is a rough estimate since the lack of symmetry in the exchange 
5k x — 8k x visible in the top panel of Fig. [T5] has not been taken into account. 

The variation of A£ with R that we obtain here is not observed in the the plane Couette flow experiments but remains 
compatible with the trend seen in CCF case [3] a]. Rather than to the role of rotation or curvature, this observation points to 
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Figure 15: Order parameter squared M 2 as a function of k x (left) and of (k x — k%) 2 (right) for L z — 48 and R = 315. 




Figure 16: Colour plot of u 2 , from left to right: R = 280, R = 300 and R = 330. L x x L z = 128 x 64. 



the role of streamwise periodic boundary conditions enforced by the cylindrical geometry or the numerical implementation. 



3.4 Dependence on R 

Variations of E, F, E%, and M against R are studied using a different protocol. Two sizes are considered: L x xL z = 128 x 64 
and 110 x 32. From the study in previous sections, both domains are expected to fit one elementary band Aj x A 2 . The 
pattern should feel "at ease" in the first domain and more "spanwise-confined" in the second one. A first simulation 
at Rq = 315 serves to prepare initial conditions for simulations at higher and lower Reynolds numbers by increasing or 
decreasing R by steps AR = 5. The flow is integrated over 5000 time units at each value of R and the so-obtained state 
is used as an initial condition for the next value of R in the range [260,350]. Additional values R = 333, 336, 337, and 
R = 370 and 390 outside the interval are also considered. At given R statistical results involve time integration over at 
least 1.5 x 10 4 time units. 

The main effect of increasing R seems to be an expansion of the turbulent part of the band pattern as illustrated in 
Fig. [16] When R is close enough to i?t the orientation of the pattern fluctuates: destroying a well-established ideal pattern, 
turbulence invades the laminar band, stays featureless for a while, before another pattern grows, which may or may not 
have the same orientation. Figure IPTl displays a featureless turbulent episode for L x x L z = 110 x 32 and R = 335, during 
which m(t) stays close to (left panel), the turbulent fraction approaches one, indicating the decrease of the size of the 
laminar domain and commanding the variation of the total energy (central panel), while the intensity of turbulence inside 
the turbulent domain does not changes (right panel). Such events cannot be mistaken with the transient occurrence of a 
defect in the pattern since both m+i and m_i remain simultaneously close to zero for a relatively long period of time, 
which is characteristic of the featureless state. They are not observed for R < 333, and go from extremely rare at R = 335 
and 337 to common at R = 340, to the most common state at R = 345 (though a trace of modulation persists). 

Since it is a three-state jump process instead of a two-state one, this feature should be treated appropriately following 



15 



E 



0.5 



1.5 




Figure 17: Time series of m (left), e and //10 (center) and et (right) zoomed on the appearance and disappearance on an 
uniformly turbulent flow for L x x L z = 110 x 32 and R = 335. 

the same procedure as for orientation fluctuations. However, it cannot be accounted for by the plain model ([3]) since 
empirical PDFs for R = 337 (Fig. [T2l right) and higher clearly present three maxima, one of which is close to the 
origin (m + ~ ~ m_). The phenomenon can however be treated within the same conceptual framework by assuming a 
slightly modified potential with an additional relative minimum at the origin separated by saddles from the main minima 
corresponding to the pattern installed in one or the other orientation, justified by the appearance of a third maximum 
in the PDFs. The splitting probability between the featureless regime and the pattern would then be controlled by the 
relative depths of the three wells [3D], which could be studied by following the procedure for orientation fluctuations [ID] . 
This complication has however not been explored further because the phenomenon is likely a size effect: In the upper 
transitional regime at large aspect-ratio, bands form out of scattered elongated regions where turbulence is depleted, sec 
Fig. 1181 The computational domains considered here arc just sufficient to contain a pattern cell of size X x x A 2 . It is 




Figure 18: Snapshot of the numerical solution for R = 340 at t = 20000; same simulation conditions as in Fig. □] (bottom- 
right), in particular L x = 432, L z — 256. 

therefore not surprising that the spatiotemporal intermittence of laminar troughs comparable in size to that cell be turned 
into temporal intermittence of well-formed laminar bands recurrently destroyed by featureless turbulence. The improved 
modelling suggested above would transform the supercritical bifurcation into a slightly subcritical one, with associated 
coexistence of featureless and patterned states, as expected from system where a spatial and temporal cohabitation of 
different states is possible. This would explain the shape of the PDFs once noise is introduced as for the original model. 
The same explanation, if correct, would explain the presence of the 'intermittent regime' described, although not fully 
investigated, by Barkley & Tuckerman [51110] since turbulence modulations around Rt are also much longer than the width 
of the oblique computational domain they considered. 

Figure [19] displays the variations of the different observables of interest with R. The growth of the width of the turbulent 
domain illustrated in Fig. [17] is clearly reflected by the increase of F with R in Fig. [191(b). Quantity F varies roughly linearly 
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Figure 19: Variation with R of the different averaged quantities for domains L x x L z — 110 x 32 (blue on line) and 128 x 64 
(red on line): (a) Perturbation energy E (full line) and turbulent energy E t (dashed line), (b) Turbulent fraction F. (c) 
Order parameter M and corresponding values from the model (dash-dotted line, <?i = 55, <?2 = 250, a — 0.002, R t — 355). 
(d) Square of the order parameter M 2 , observations (o : 128 x 64 and + : 110 x 32) and linear fits (lines). 



with R for R < R t , and with a much smaller slope above. This increase mostly explains the growth of the perturbation 
energy since the turbulent energy depends more weakly on R, with no singular behaviour visible at R t (Fig. 1191 a): the 
intensity of turbulence inside the turbulent domains does not seem sensitive to the global organisation in an oblique pattern. 
The slope discontinuity at Rt — 350 marking the bifurcation was used as a criterion in our previous study [13] • Values 
obtained for E, E t , and F are slightly different for the two sizes considered, which is related to lateral confinement effects 
already illustrated in Fig. [13] (left) for R = 315. 

In Figure fT9l (c) . it can be seen the order parameter M departs from the expected classical e 1 / 2 behaviour and tends to 
saturate in the lowest part of the transitional range. For L x x L z = 110 x 32, it even decreases as R is lowered further, 
which is again a confinement effect since, from the experiments [3] as well as from our earlier (less well resolved) numerical 
results [13], the spanwise wavelength X z is expected to increase up to about 80 as R decreases: this implies a less optimal 
pattern and a weaker modulation for L z = 32, while for L x x L z = 128 x 64, with a more favourable A z , M continues to 
increase as R is lowered in agreement with the Ginzburg-Landau picture. In the upper part of the transitional range, M 
decreases quickly as R increases. The decay of the modulation corresponds to the increase of the width of the turbulent 
domain. Again in line with the Ginzburg-Landau interpretation, the variation of M 2 with R (Fig. [121 d) appears to be 
linear with a slope l/(g\Rt) — —5.3 x 10 -5 . Meanwhile, the extrapolation of M 2 to zero gives a value of Rt ~ 355 or 348 
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depending on whether one takes the data from case L x x L z = 110 x 32 or 128 x 64, respectively. In contrast with what 
happens for R ~ R g , here the estimate with L z = 32 is likely the best one since experiments suggest X z ~ 37.5 either by 
extrapolation for plane Couette flow or from measurements in the CCF [3 . Taking R t ~ 355, we get g\ ~ 55, which is 
consistent with Prigent's value g\ ~ 100 in the CCF case. This value of g\ further yields £ z ~ 2.3 (measured values for 
CCF range between 0.9 and 3.2) and ( x ~ 11 for R = 315. 

In Fig. [19] (c), it can be noticed that M remains finite for R > Rt, as the result of intrinsic fluctuations in the 
featureless turbulent regime, in contrast with what would happen in the deterministic case. Fluctuations indeed gives a 
finite background level to modes m±, a fact which is well accounted for by the model in the mean- field approximation 
represented by dash-dotted line in Fig [IS] (c). 

These result are not qualitatively affected by the increase of resolution from N y = 15 and N x , z /L XtZ = 8/3 to N y = 27 
and N XtZ /L x , z — 4 as can be seen in Figure [201 which compares the results for both resolutions at size L x x L z — 110 x 32. 
The quantitative change is minor and in the expected fashion [13J . The thresholds R g and Rt move to approximately 321 
and 390. The square of the perturbation undergoes an increase of about 10%. Apprt from the threshold shift, the turbulent 
fraction F is little affected by the resolution change. Both E and F display the expected slope-break. Given the uncertainty 
on the values of M near R t , the value gi ~ 30 at N y — 27 is acceptable. This reassert the validity of our semi-quantitative 
approach. 

Parameter g-i has little influence and reasonable results are obtained from 0.001 < a < 0.003. This estimate is consistent 
with the value obtained from the fit of the phase dynamics fit afro ~ 4 x 10~ 4 , (Fig. E} if we accept Prigent's finding 
To ~ 30h/U p2 (a,c)]. The variance of the fluctuations of m in the vicinity of Rt is also of interest. Let us define: 

S 2 = 2 I {m! - M) 2 n(m', m") Am! dm" . 

J m" <m ; 

Figure [21] displays the variation of S as a function of R. Fluctuations appear to be strongly enhanced in the vicinity of 
Rt, due to orientation changes and re-entrance of featureless turbulence. Though model ([3]) does not account for the latter 
phenomenon, it already explains a large part of the enhancement. Its parameter gi controls the amplitude of fluctuations 
that bring about orientation changes. The position of the maximum of S(R) strongly depends on it. With g\ = 55, 
satisfactory agreement is found for gi > 120. Results obtained with 52 = 250 are represented as a dash-dotted line in 
Fig. 1211 Including the re-entrance of featureless turbulence would certainly increase the variability but this would still not 
be the whole story since, like for second order phase transitions, one would expect a divergence of S in the form S oc |e| 7 , 
7 being the critical exponent attached to the susceptibility of the order parameter, just rounded off by finite-size effects. 
Even at reduced numerical resolution, improving the statistics to study the pattern's fluctuations in the simulations seems 
presently out of reach. 
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Figure 21: Average root-mean-square fluctuation of m as a function of R for L x x L z = 110 x 32 (full line), 128 x 64 
(dashed), and the model with gi = 55, gi = 250, R t = 355, a = 0.002 (dash-dotted). 

4 Summary and Conclusion 

Prigent et al. [5] have put the problem of the emergence of turbulent bands in wall-bounded flows within the Ginzburg- 
Landau framework of pattern formation, adding noise to account for background turbulence. Doing so, they were able to 
extract most of the coefficients in the model equation from laboratory experiments in the case of circular Couette flow, 
while restricting themselves to threshold localisation and wavelength measurements for PCF. In a similar vein, Barkley et 
al. [10] later performed simulations of PCF, detecting the formation of bands from Fourier analysis of the pattern. They 
considered a quasi-one-dimensional configuration excluding orientation fluctuations expected to play a role close to Rt for 
symmetry reasons. Though having the model in mind, they did not attempt any quantitative fit. Our work has been 
mostly devoted to overcome these two limitations, to check the validity of the noisy Ginzburg-Landau framework, and to 
compare finding for PCF to those for CCF. Previous result [T3] were reasserted, showing that controlled under-resolution 
gives excellent qualitative agreement with experiments and good quantitative results once corrected for a general shift of 
the range [R g , R t ] where the bands are present. We performed numerical experiments in domains of sizes able to contain 
one to three bands in the spanwise direction and one or two bands in the streamwise direction, while letting the pattern's 
orientation fluctuate. Under-resolution reducing the computational load, we could carry out long duration simulations in 
order to accumulate reliable statistics. 

The emergence of bands was first quantitatively characterised using standard statistical quantities such as the total 
perturbation energy E, the turbulent fraction F, and the average energy contained in turbulent domains E t . These 
quantities quickly converge to their steady-state values but do not give information on orientation or wavelength fluctuations. 
This limitation has been next overcome by defining order parameters measuring the amplitude of the modes involved in the 
Fourier series decomposition of the patterns, appropriately amending the Barkley et al. definitions and procedure. The full 
nonlinear dispersion relation describing the formation of bands could be studied by varying the Reynolds number and the 
size of the computational domain which controls the allowed wavevectors. The coefficients of the relevant Ginzburg-Landau 
equation and the intensity of the noise were estimated, showing the overall consistency of the approach. In particular, two 
coherence lengths, spanwise and streamwise, were evaluated and the square of the modulation amplitude was shown to vary 
linearly with R far enough from R t , while its fluctuations and the intermittent re-entrancc of featureless turbulence were 
strongly enhanced close to R t . It has been argued that the re-entrance of featureless turbulence was a side effect of the 
limited size of the system, probably explaining the 'intermittent regime' of Barkley & Tuckerman [5] by the same token, 
and that this observation should be better replaced on a spatiotemporal footing in more extended domain, in relation 
to patterns with mixed orientations observed near Rt in CCF experiments 3^ or PCF simulations in Fig. [18] Finally, 
comparing our results with those obtained in CCF we obtain satisfactory general agreement, but with the supplementary 
information that the streamwise coherence length £ x is significantly larger than the spanwise coherence length £ z indicating 
that the selection of the streamwise wavelength X x is more effective than that of the spanwise wavelength A z . 
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As a whole, the emergence of oblique bands from featureless turbulence upon decreasing R has been seen to fit the 
conventional framework of a pattern-forming instability. However, the very fact that the base state is turbulent calls for the 
introduction of a large noise in the picture. These numerical studies are performed with the hope that they will contribute 
to the understanding of the cohabitation of turbulent and laminar flow typical of the transition to/from turbulence in 
wall-bounded flows, the detailed mechanism of which is still largely unknown. 
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